Blowing up Feynman integrals 



Christian Bogner a and Stefan WeinzierP 

a Institut fur Physik, Universitat Mainz, D - 55099 Mainz, Germany 

In this talk we discuss sector decomposition. This is a method to disentangle overlapping singularities through 
a sequence of blow-ups. We report on an open-source implementation of this algorithm to compute numerically 
the Laurent expansion of divergent multi-loop integrals. We also show how this method can be used to prove 
a theorem which relates the coefficients of the Laurent series of dimensionally regulated multi-loop integrals to 
periods. 



1. INTRODUCTION 

The calculation of higher-order corrections in 
perturbation theory in particle physics relies to 
a large extent on our abilities to compute Feyn- 
man loop integrals. These calculations are com- 
plicated by the occurrence of ultraviolet and in- 
frared singularities. Ultraviolet divergences are 
related to the high-energy behaviour of the inte- 
grand. Infrared divergences may occur if massless 
particles are present in the theory and are related 
to the low-energy or collinear behaviour of the in- 
tegrand. 

Dimensional regularisation is usually employed 
to regularise these singularities. Within dimen- 
sional regularisation one considers the loop inte- 
gral in D space-time dimensions instead of the 
usual four space-time dimensions. The result 
is expanded as a Laurent series in the parame- 
ter e — (4 — D)/2, describing the deviation of 
the D-dimensional space from the usual four- 
dimensional space. The singularities manifest 
themselves as poles in 1/e. Each loop can con- 
tribute a factor 1/e from the ultraviolet diver- 
gence and a factor 1/e 2 from the infrared diver- 
gences. Therefore an integral corresponding to a 
graph with I loops can have poles up to 1 /e 21 . 

At the end of the day, all poles disappear: The 
poles related to ultraviolet divergences are ab- 
sorbed into renormalisation constants. The poles 
related to infrared divergences cancel in the final 
result for infrared-safe observables, when summed 
over all degenerate states or are absorbed into 



universal parton distribution functions. The sum 
over all degenerate states involves a sum over con- 
tributions with different loop numbers and differ- 
ent numbers of external legs. 

However, intermediate results are in general a 
Laurent series in e and the task is to determine 
the coefficients of this Laurent series up to a cer- 
tain order. Sector decomposition [1-6] can be 
used to compute numerically the coefficients of 
the Laurent series. At the same time, sector de- 
composition provides the tools to prove that un- 
der rather weak assumptions the coefficients of 
the Laurent series are special numbers called pe- 
riods [7]. 

2. MULTI-LOOP INTEGRALS 

Let us now consider a generic scalar Z-loop in- 
tegral Iq in D = 2m — 2e dimensions with n 
propagators, corresponding to a graph G. Let 
us further make a slight generalisation: For each 
internal line j the corresponding propagator in 
the integrand can be raised to a power Vj . There- 
fore the integral will depend also on the numbers 
z/ lr ..,z/ n . It is sufficient to consider only the case, 
where all exponents are natural numbers: Vj £ N. 
We define the Feynman integral by 
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with v = V\ + ... + v n . The momenta qj of the 
propagators are linear combinations of the exter- 
nal momenta and the loop momenta. The pref- 
actors are chosen such that after Feynman pa- 
rameterisation the Feynman integral has a simple 
form: 
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The functions U and J 7 depend on the Feynman 
parameters and can be derived from the topology 
of the corresponding Feynman graph G. Cutting I 
lines of a given connected Hoop graph such that it 
becomes a connected tree graph T defines a chord 
C(T, G) as being the set of lines not belonging to 
this tree. The Feynman parameters associated 
with each chord define a monomial of degree I. 
The set of all such trees (or 1-trees) is denoted 
by T\. The 1-trees T G Ti define U as being 
the sum over all monomials corresponding to the 
chords C(T,G). Cutting one more line of a 1- 
tree leads to two disconnected trees (Ti,T-z), or 
a 2-tree. T2 is the set of all such pairs. The 
corresponding chords define monomials of degree 
I + 1. Each 2-tree of a graph corresponds to a cut 
defined by cutting the lines which connected the 
two now disconnected trees in the original graph. 
The square of the sum of momenta through the 
cut lines of one of the two disconnected trees T\ 
or T2 defines a Lorentz invariant 
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The function To is the sum over all such monomi- 
als times minus the corresponding invariant. The 
function T is then given by To plus an additional 
piece involving the internal masses nij. In sum- 
mary, the functions U and T are obtained from 



the graph as follows: 
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In eq. ([T]) we just considered scalar integrals. A 
priori more complicated cases, where the loop mo- 
mentum appears in the numerator might occur. 
However, there is a general reduction algorithm, 
which reduces these tensor integrals to scalar in- 
tegrals [8,9]. The price we have to pay is that 
these scalar integrals involve higher powers of the 
propagators and/or have dimensions shifted by 
two units. Therefore we introduced in eq. ([1]) 
higher powers of the propagators and kept the 
dimension D = 2m — 2e arbitrary. As a conse- 
quence, integrals of the form as in eq. {1} are the 
most general loop integrals we have to consider. 

3. SECTOR DECOMPOSITION 

In this section we review the algorithm for it- 
erated sector decomposition. The starting point 
is an integral of the form 

/n / n \ r 

(5) 

where fa = a,; + ebi and Xj = cj + edj . The inte- 
gration is over the standard simplex. The a's, 6's, 
c's and d's are integers. The P's are polynomials 
in the variables x%, x n . The polynomials are 
required to be non-zero inside the integration 
region, but may vanish on the boundaries of the 
integration region. The algorithm consists of the 
following steps: 

Step 0: Convert all polynomials to homogeneous 
polynomials. 

Step 1: Decompose the integral into n primary 
sectors. 
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name "blow-up" . 

Step 3: The singular behaviour of the integral 
depends now only on the factor 



(9) 



We Taylor expand in the integration variables and 
perform the trivial integrations 



Figure 1. Illustration of sector decomposition and 
blow-up for a simple example. 
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leading to the explicit poles in 1/e. 



Step 2: Decompose the sectors iteratively into 
sub-sectors until each of the polynomials is of the 
form 



P 
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where c ^ and P' '(x) is a polynomial in the vari- 
ables Xj without a constant term. In this case the 
monomial prefactor x™ 1 ...x™ n can be factored 
out and the remainder contains a non-zero con- 
stant term. To convert P into the form ([6|) one 
chooses a subset S = {a\, ak} C {1, ...n} 
according to a strategy discussed in the next sec- 
tion. One decomposes the fc-dimensional hyper- 
cube into k sub-sectors according to 

l d n x = J2 [ II 6 ( x <*i > • (?) 
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In the l-ih sub-sector one makes for each element 
of S the substitution 



for i ^ I. 



(8) 



This procedure is iterated, until all polynomials 
are of the form (J6j) . Fig. Q] illustrates this for 
the simple example S = {1,2}. Eq. Q gives the 
decomposition into the two sectors x\ > x% and 
X2 > x±. Eq. |8]) transforms the triangles into 
squares. This transformation is one-to-one for all 
points except the origin. The origin is replaced 
by the line X\ = in the first sector and by the 
line X2 = in the second sector. Therefore the 



Step 4: All remaining integrals are now by con- 
struction finite. We can now expand all expres- 
sions in a Laurent series in e and truncate to the 
desired order. 

Step 5: It remains to compute the coefficients 
of the Laurent series. These coefficients contain 
finite integrals, which can be evaluated numeri- 
cally by Monte Carlo integration. 

4. HIRONAKA'S POLYHEDRA GAME 

In step 2 of the algorithm we have an iteration. 
It is important to show that this iteration termi- 
nates and does not lead to an infinite loop. There 
are strategies for choosing the sub-sectors, which 
guarantee termination. These strategies [10-14] 
are closely related to Hironaka's polyhedra game. 

Hironaka's polyhedra game is played by two 
players, A and B. They are given a finite set M 
of points m = (mi, m n ) in N™ , the first quad- 
rant of N". We denote by A c WL the positive 
convex hull of the set M. It is given by the convex 
hull of the set 
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The two players compete in the following game: 

1. Player A chooses a non-empty subset S C 
{1, ...,«}. 

2. Player B chooses one element i out of this 
subset S. 
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Figure 2. Illustration of Hironaka's polyhedra game. 

Then, according to the choices of the players, the 
components of all (mi, m n ) € M are replaced 
by new points (m^, m' n ), given by: 

m'j = nij, if j^i, 

where for the moment we set c = 1. This defines 
the set M'. One then sets M — M' and goes 
back to step 1. Player A wins the game if, after 
a finite number of moves, the polyhedron A is of 
the form 

A = to + M™, (13) 

i.e. generated by one point. If this never occurs, 
player B has won. The challenge of the polyhedra 
game is to show that player A always has a win- 
ning strategy, no matter how player B chooses his 
moves. A simple illustration of Hironaka's poly- 
hedra game in two dimensions is given in fig. [U 
Player A always chooses S = {1, 2}. In ref. [5] we 
have shown that a winning strategy for Hiron- 
aka's polyhedra game translates directly into a 
strategy for choosing the sub-sectors which guar- 
antees termination. 

5. AN OPEN-SOURCE IMPLEMENTA- 
TION 

The algorithm for computing numerically the 
coefficients of the Laurent series of a divergent 
multi-loop integral involves symbolic computer 
algebra and numerical integration. Partially due 
to this "double nature" of the algorithm, no pub- 
lic program was available until recently. We im- 



plemented3 the algorithm in the language C++ 
using the library GiNaC [15] for the symbolic 
part. In addition the GNU Scientific Library 
[16] is used for the numerical Monte Carlo in- 
tegration. We have implemented various strate- 
gies for choosing the sub-sectors. Three of them 
are guaranteed to terminate. However in prac- 
tice they lead to a rather large number of sub- 
sectors. In addition we have implemented two 
heuristic strategies, for which we neither have 
a proof that they terminate, nor do we know a 
counter-example which would lead to an infinite 
recursion. If these heuristic strategies terminate, 
the number of the generated sub-sectors tends to 
be smaller than the corresponding numbers for 
the other strategies. 

6. PERIODS 

Periods are special numbers. Before we give the 
definition, let us start with some sets of numbers: 
The natural numbers N, the integer numbers Z, 
the rational numbers Q, the real numbers R and 
the complex numbers C are all well-known. More 
refined is already the set of algebraic numbers, de- 
noted by Q. An algebraic number is a solution of 
a polynomial equation with rational coefficients: 

x n + a n - lX n - 1 + • • • + a = 0, G Q. (14) 

As all such solutions lie in C, the set of algebraic 
numbers Q is a sub-set of the complex numbers C. 
Numbers which are not algebraic are called tran- 
scendental. The sets N, Z, Q and Q are countable, 

1 The program can be obtained from 
http : //www . higgs . de/~ st ef anw/sof twar e . html 
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whereas the sets R, C and the set of transcenden- 
tal numbers are uncountable. 

Periods are a countable set of numbers, lying 
between Q and C. There are several equivalent 
definitions for periods. Kontsevich and Zagier 
gave the following definition [17]: A period is a 
complex number whose real and imaginary parts 
are values of absolutely convergent integrals of 
rational functions with rational coefficients, over 
domains in 1" given by polynomial inequalities 
with rational coefficients. Domains defined by 
polynomial inequalities with rational coefficients 
are called semi-algebraic sets. 

We denote the set of periods by P. The alge- 
braic numbers are contained in the set of periods: 
Q £ P. In addition, P contains transcendental 
numbers, an example for such a number is ir: 



dx dy. 



(15) 



x 2 -\-y 2 <l 



The integral on the r.h.s. clearly shows that ir is a 
period. On the other hand, it is conjectured that 
the basis of the natural logarithm e and Euler's 
constant je are not periods. Although there are 
uncountably many numbers, which are not peri- 
ods, only very recently an example for a number 
which is not a period has been found [18]. 

We need a few basic properties of periods: The 
set of periods P is a Q-algebra [17,19]. In partic- 
ular the sum and the product of two periods are 
again periods. 

The defining integrals of periods have inte- 
grands, which are rational functions with ratio- 
nal coefficients. For our purposes this is too re- 
strictive, as we will encounter logarithms as inte- 
grands as well. However any logarithm of a ra- 
tional function with rational coefficients can be 
written as 



lng(x) 



dt 



(g(x)-l)t + l 



(16) 



7. A THEOREM ON FEYNMAN INTE- 
GRALS 

Let us consider a general scalar multi-loop in- 
tegral as in eq. @ Let m be an integer and set 



D = 2m — 2e. Then this integral has a Laurent 
series expansion in e 



j=-2l 

Theorem 1: In the case where 



(17) 



1. all kinematical invariants st are zero or neg- 
ative, 

2. all masses rrii and \i are zero or positive 
(,i^0), 

3. all ratios of invariants and masses are ratio- 
nal, 

the coefficients Cj of the Laurent expansion are 
periods. 

In the special case were 

1. the graph has no external lines or all invari- 
ants st are zero, 

2. all internal masses nij are equal to fi, 

3. all propagators occur with power 1, i.e. 
Vj = 1 for all j, 

the Feynman parameter integral reduces to 
„ n 
I G = d n xS(l-Y,^)U~ D/2 (18) 

Xj >o 1=1 

and only the polynomial U occurs in the inte- 
grand. In this case it has been shown by Belkale 
and Brosnan [20] that the coefficients of the Lau- 
rent expansion are periods. 

Using the method of sector decomposition we 
are able to prove the general case. We will actu- 
ally prove a stronger version of theorem 1. Con- 
sider the following integral 
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The integration is over the standard simplex. The 
a's, 6's, d's and fs are integers. The P's are 
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polynomials in the variables X\, x n with ra- 
tional coefficients. The polynomials are required 
to be non-zero inside the integration region, but 
may vanish on the boundaries of the integration 
region. To fix the sign, let us agree that all poly- 
nomials are positive inside the integration region. 
The integral J has a Laurent expansion 

oo 

J = ]Tc^. (20) 

3=30 

Theorem 2: The coefficients Cj of the Laurent 
expansion of the integral J are periods. 

Theorem 1 follows then from theorem 2 as the 
special case a* = Vi — 1, 6j = 0, r = 2, Pi = U, 
P 2 = T, d x + eh = v-{l + l)D/2 and d 2 + ef 2 = 
ID/2 - v. 

Proof of theorem 2: To prove the theorem we 
will give an algorithm which expresses each co- 
efficient Cj as a sum of absolutely convergent in- 
tegrals over the unit hypercube with integrands, 
which are linear combinations of products of ra- 
tional functions with logarithms of rational func- 
tions, all of them with rational coefficients. Let 
us denote this set of functions to which the inte- 
grands belong by M.. 

The unit hypercube is clearly a semi-algebraic 
set. It is clear that absolutely convergent in- 
tegrals over semi-algebraic sets with integrands 
from the set M are periods. In addition, the sum 
of periods is again a period. Therefore it is suffi- 
cient to express each coefficient Cj as a finite sum 
of absolutely convergent integrals over the unit 
hypercube with integrands from M. To do so, 
we use iterated sector decomposition. Details of 
the proof can be found in [7] . 

8. CONCLUSIONS 

In this talk we reported on an open-source im- 
plementation of a computer program to compute 
numerically the coefficients of the Laurent expan- 
sion of divergent multi-loop integrals. We also 
showed that under rather weak assumptions all 
these coefficients are periods. The technique we 
used was based on sector decomposition and se- 
quences of blow-ups. It appears that Feynman in- 
tegrals and algebraic geometry are not un-related, 
as was already pointed out in [21]. 
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